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We review the calculation of moments of both the polarized and unpolarized parton 
distribution functions of the nucleon in lattice QCD, and in particular their extrapolation 
to the physical region. We also discuss the reconstruction of the x dependence of the 
valence quark distributions in the nucleon from a finite number of lattice moments. 

1. Introduction 

One of the defining features of any hadron are its quark and gluon (or generically, 
parton) momentum distributions. Considerable information has been accumulated 
on parton distribution functions (PDFs) of the nucleon and nuclei from deep inelas- 
tic scattering and other high energy experiments. Information from the Drell-Yan 
reaction in irN scattering has also been used to determine PDFs of the pion. 

Theoretically, the operator product expansion (OPE) in QCD allows one to 
isolate the soft, non-perturbative PDFs from the hard, perturbatively calculable 
scattering processes. This factorization property gives rise to a universal set of 
PDFs which can be used to describe various reactions. Because PDFs parameterize 
the interactions of quarks and gluons in the hadron over all distance scales, they 
are sensitive to the dynamics of QCD in the strong coupling regime - i.e., at scales 
Q 2 ~ Aq CD . The only rigorous theoretical approach currently capable of calculating 
observables in the non-perturbative regime from first principles is lattice QCD. 

The first exploratory studies of PDFs on the lattice were made almost two 
decades agoP Advances in computing power have recently enabled dedicated sim- 
ulations to be undertaken which for the first time can be directly compared with 
phenomenologyPEEEl Since the numerical simulations are performed in Euclidean 
space-time, it is not possible to compute the PDFs directly as a function of the 
light-cone momentum fraction, x (PDFs are formally defined as light-cone corre- 
lation functions involving currents with space-time separation z 2 — (ct) 2 ~ 0). On 
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the other hand, using the OPE one can formally express the moments of PDFs in 
terms of matrix elements of local operators between hadron states. These matrix 
elements can be calculated numerically on the lattice. 

Despite the impressive recent progress in the numerical simulations, several ap- 
proximations are still necessary in order to relate the lattice moments to experiment. 
Firstly, since space-time is discretized on the lattice, with some finite lattice spac- 
ing a, the results must be extrapolated to the continuum limit, a — > 0. The finite 
number of lattice sites also means that an extrapolation to the infinite volume limit 
is necessary to avoid omitting important physics arising from the long-range part of 
the nucleon wave function. In addition, since the computational cost of simulations 
scales with the quark mass roughly as m~ , current lattice simulations are performed 
with unphysically large masses for u and d quarks, typically m^ att > 30 MeV, so 
that an extrapolation to physical masses, m^ hys w 5 MeV, is essential.^ 

The chiral (m q ) extrapolation is found to play a particularly vital role in under- 
standing the connection between the lattice results and phenomenology. Whereas 
the lattice simulations yield results for the moments of unpolarized quark dis- 
tributions in the nucleon which are typically 50% larger than experiment when 
extrapolated linearly to mj tys , inclusion of the non-linear dependence on m q aris- 
ing from the long-range structure of the nucleon removes most of the discrepancy.^ 
Furthermore, significant finite volume effects have been found"! in the calculation of 
the axial vector charge gA, which may eventually explain the residual discrepancy 
between the extrapolated lattice value and experiment. 

In this mini-review we focus on the challenges involved with extrapolating the 
lattice data from the currently accessible regions of parameter space to the physical 
region. Using constraints available from chiral effective theory, in Sec. 2 we review 
extrapolation formulas which seek to describe the lattice data over a large range of 
quark masses, and in the chiral limit. The moments of both the spin-averaged and 
spin-dependent quark distributions in the nucleon are considered (a corresponding 
analysis of the moments of the quark distributions in the pion was also undertaken 
recently 21 ) . The more ambitious task of reconstructing the x dependence of PDFs 
from a finite number of moments is discussed in Sec. 3. Because of the small number 
of lattice moments available, only valence distributions can be analyzed at present; 
sea quark distributions await future lattice data on higher moments of PDFs which 
are necessary for an independent determination of the valence and sea components. 
Finally, in Sec. 4 we draw some conclusions, and outline prospects for the recon- 
struction of PDFs from future lattice simulations. 
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2. Extracting physical results from lattice simulations 

The moments of parton distributions in the nucleon are formally defined as 



for the spin-averaged, helicity and transversity distributions, respectively. From 
their definition, the moments alternate between the total (q + q) and valence (q — q) 
distributions, depending on whether n is even or odd. Using the OPE, these mo- 
ments can be related to ground state hadron matrix elements of specific twist-two 
operators,^ which are calculated on the lattice. 

In this review we focus not on the details of the lattice simulations of the mo- 
ments (a comprehensive survey of results was recently given by the author ^M^) , but 
on their physical interpretation. In particular, we examine the consequences of chi- 
ral symmetry for the extrapolation of PDF moments to physical values of the quark 
mass, m q . 

2.1. Chiral symmetry 

The importance of c hiral sy mmetry and the role of the pion cloud in hadronic 
fTTTT^T 

physics is well known! 11 At small quark masses, hadronic observables can be 
systematically expanded in a series in m q within the framework of chiral perturba- 
tion theory (%PT). While the expansion coefficients are generally free parameters, 
one of the unique consequences of spontaneous chiral symmetry breaking in QCD is 
the appearance of terms involving odd powers or logarithms of m, . From the Gell- 
Mann-Oakes-Renner relation, which relates the quark and pion masses at small 
rn-n, ~ m q , one finds that such terms are non-analytic in the quark mass. Fur- 
thermore, their coefficients, which are determined from the infrared behavior of the 
pion loops, are generally model independent. 

The non-analytic term involving the lowest power of is known as the "leading 
non-analytic" (LNA) term in the expansion. For the moments of the nucleon PDFs 
this was shown by Thomas et aZP*21 to have the generic behavior m\ log arising 
from ttN intermediate states. This was later confirmed in yPT, w here the coefficients 
of these terms were also calculated, both for the nucleonp^l and for the pionP^ 
Using these constraints, a low order chiral expansion for the moments of the PDFs 
was developed recently by Detmold et al. I 7 I 8 I 1U I 16 I 17 I which incorporated the LNA 
behavior of the moments as a function of m q and also provided a connection to the 
heavy quark limit. 




(1) 
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2.2. Chiral extrapolation of PDF moments 

Direct comparison of the currently available lattice moments of nucleon PDFs with 
phenomcnological ones is non-trivial. The lattice moments in general receive con- 
tributions from diagrams in which the local operator insertions are on quark lines 
which arc connected to the nucleon source, as well as those where the operator acts 
on a quark loop in the vacuum, which is connected only by gluon lines to the quark 
lines originating in the nucleon source. Obtaining a signal from the latter, so- call ed 
"disconnected graphs" is extremely difficult on the lattice, as recent studies^! of 
the strangeness form factors of the nucleon have shown, and to date only the con- 
nected contributions to the lattice PDF moments have been computed. Fortunately, 
the near degeneracy of the u and d quark masses, which makes charge symmetry 
such a good symmetry of nature, means that the disconnected diagrams cancel 
very precisely in the isovector difference between the u and d PDFs of the proton. 
An unambiguous comparison of the connected isovector PDF lattice moments can 
therefore be made with the analogous experimental moments. 

The mass dependence of the mo ments of the unpolarized isovector PDF of the 
nucleon can be parameterized a 



1 + CLNA"4 2 r 2 + 5 " 2/2 » ( 4 ) 

where (for n > 0) the chiral coefficient is given fy |14 | 15 | CLNA _ — / (47r^7r) 2 . 
Although we shall work only with full QCD, we note that the corresponding coeffi- 
cients have al so b een computed in quenched and partially-quenched chiral pertur- 
bation theorypj The n = moment, which corresponds to isospin charge, is not 
renormalized by pion loops. The argument of the logarithm corresponds to the case 
where the pion loops are regulated with a simple, sharp three-momentum cut-off. 
The parameter [i has the effect of suppressing the rapid variation of the chiral log- 
arithm away from the chiral limit. Physically it is related to the size of the nucleon 
core, which acts as the source of the pion fieldP Since the isovector distribution 
u(x) — d(x) — > S(x — 1/3) in the heavy quark limit (i.e., as m q — > 00), one may 
choose to constrain the constant b n in Eq. by 

b n = — - a n (l - /i 2 c L na) • (5) 

The results of the fits to the lattice data for the n = 1 , 2 and 3 moments of the 
unpolarized u — d distribution using Eq. are shown in Fig. 1. In each plot the 
central curve is shown with two error bands: the inner (darkly shaded) band shows 
variation of fi by ± 20%, whilst the outer (lightly shaded) shows the additional 
effects that result from shifting the lattice data up and down within the extent of 
their error bars. The extrapolated moments at the physical pion mass are listed in 
Table [U The fits to the data are quite insensitive to the choice of m^n (as long as it 
is large) and it has been set to 5 GeV for all n. In practice, fits to the lattice data 
in which the heavy quark limit is not guaranteed (i.e., where the term involving 
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Fig. L Moments of the unpolarized u — d distribution in the proton, for n = 1, 2 and 3. Lattice 
dat £ . 10 include both quenched (solid symbols) and unquenched (open symbols) results. The solid 
line represents the full chiral extrapolation, while the inner (darkly shaded) error band shows 
variation of fj. by ± 20%, with the outer band (lightly shaded) showing the additional effects of 
shifting the lattice data within the extent of thei r er ror bars. Linear extrapolations are indicated 
by dashed lines, and the phenomenological valueJ^Hl are shown as large stars at the physical pion 
mass. 



b n is simply 6„mJ) and b n is a third fitting parameter,^ are indistinguishable from 
those in Fig. 1. 

Note that the majority of the data points (filled symbols) are obtained from 
simulations employing the quenched approximation (in which background quark 
loops are neglected) whereas Eq. (@J is based on full QCD with quark loop effects 
included. On the other hand, recent calculations with dynamical quarks suggest that 
at the relatively large pion masses (m w > 0.5-0.6 GeV) where the full simulations 
are currently performed, the effects of quark loops are largely suppressed, as the data 
in Fig. 1 (small open symbols) indicate. Further details of the lattice datap^EEl 
and a more extensive discussion of the fit parameters, can be found elsewhere!^ 

A similar analysis leads to analogous lowest o rder LNA parameterizations of the 
mass dependence of the spin-dependent moments^! 



and 



)Au-Ad = Aq„ 1 + AcLNAffl^ log ~ — p — r ) + Ab n - 2 , (6) 



{x n )su-Sd = Sa n [ 1 + Sc LNA mi log 2 * 2 ) + Sb n _ 2 ; ^ 2 , (7) 
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Tabic 1. Experime ntal momentsSOEH and 
moments extrapolatecED from lattice data. 



Moment 


Experiment 


Extrapolated 


(x L )u-d 


0.145(4) 


0.18(3) 


{X 2 )u-d 


0.054(1) 


0.05(2) 


{x Z )u-d 


0.022(1) 


0.02(1) 


(l°)A«-Ad 


1.267(4) 


1.12(6) 


{x 1 )Au-Ad 


0.210(25) 


0.27(2) 


{x 2 )Au-Ad 


0.070(11) 


0.14(4) 


(l)su-Sd 




1.23(7) 


{x)Su-Sd 




0.51(9) 



where the LNA coefficients are given by^' Aclna = — (1 + 2g^)/(47r/ w ) 2 and 
^clna = — (1 + 4.9a) / ^2(47r -J f 7r ) 2 . In the heavy quark limit, where spin -flavor sym- 



metry is exact, both Au(x) — Ad(x) and 5u(x) — 5d(x) are given b 
which leads to the constraints 



Ab„ 



3 rH 



and 



3 ™+i 



Aa n (l - [i 2 Ac LNA ) , 



Sa n (l - ^ 2 5c LNA ) 



\S(x-l/3), 



(8) 



(9) 



These are the most general lowest order parameterizations of the twist- 2 PDF mo- 
ments consistent with chiral symmetry and the heavy quark limits of QCD. 

The n = moment of the spin-dependent isoyector PDF is equivalent to the 
axial charge of the nucleon, gA- It is well knownpSfor instance through the Adler- 
Weisberger sum rule, that the A resonance plays a key role in gA- In the framework 
of the chiral expansion, vertex renormalization contributions involving a A isobar, 
although not leading non-analytic, are lar ge and extremely important in countering 
the effect of wave function renormalization! Indeed , it is entirely because of the 
explicit role of the A that the bare and renormalized irNN couplings are closePS 

From the point of view of the chiral extrapolation problem, the tendency for 
cancellation between wave function and polarized vertex corrections means that, 
once virtual it A loops are included, the chiral variation as — > is approximately 
linear. That is, for the polarized moments the chiral non-analytic behavior produces 
far less curvature as one approaches the physical pion mass than in the unpolarized 
case. This can be seen explicitly in Figs. 2 and 3, where the n = 0, 1 and 2 moments 
of the helicity distribution Au(x) — Ad(x) and the n = and 1 moments of the 
transversity distribution Su(x) — Sd(x) are plotted, respectively. In these figures, the 
long-dashed curves show the extrapolations ignoring the irA contributions. 

The curvature in the extrapolated lowest moments of the helicity distribution 
(gA) increases somewhat the discrepancy between the lattice value and experiment, 
although the contributions of the A largely reduce this effect. Nevertheless, there 
does appear to be a residual 10-15% underestimation of gA, while at the same 
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Fig. 2. Moments of the helicity distribution Am — Ad in the nucleon for n = 0, 1 and 2. The 
long-dashed curves show extrapolations neglecting the A. The meaning of all other symbols is as 
in Fig. 1. 



time the n = 1 and 2 moments of Au — Ad are somewhat overestimated (although 
the errors on the n = 2 lattice data in particular are rather large). There have 
been several suggestions for the possible origin of this discrepancy, including the 
effects of working with a finite volume on the lattice. The RBCK Collaboration 
has in fact observed an unusually strong dependence of gA on the lattice volume!^ 
This is illustrated in the top panel of Fig. [21 where the large open boxes and dia- 
monds correspond to calculations on lattices with spatial volumes of (2.4 fm) and 
(1.2 fm) 3 , respectively. Clearly the issue of understanding the value of gA in QCD 
is of central importance in hadronic physics; moreover, resolution of this apparent 
discrepancy should enable a more reliable prediction to be made for the moments 
of the transversity distribution in Fig. 3, for which there are as yet no experimental 
data. 



2.3. Volume dependence 

The origin of the strong lattice volume dependence of gA is likely related to the 
fact that the pion held does not vanish before it meets the lattice boundary!^ 
However, there is another important issue to consider when it comes to the curvature 
that we have been discussing with regard to the chiral extrapolation problem. One 
would really like to see confirmation of this curvature as the lattice calculations are 
continued to smaller pion mass and eventually use these calculations to determine 
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Fig. 3. Moments of the isovector transversity distribution Su — Sd for n = and 1. Symbols and 
curves are as in Fig. 2. 




0.1 1 ■ ■ ■ ' ■ ■ ■ ' ■ ■ ■ ' ■ ■ ■ 1 

0.2 0.4 0.6 0.8 

n% \_GeV 2 ] 

Fig. 4. Estimate of the volume dependence of the chiral curvature in the extrapolation of the 
isovector quark momentum fraction, (x) u _a, obtained with a simple cut-off on the lower limit of 
the pion three momentum, for various box side lengths L. 

the parameter \i. This is difficult for two reasons. First, the chiral coefficients are 
smaller in quenched simulations (which are numerically less demanding), and so the 
curvature may not be observable until very low masses are reached. Secondly, even 
in the case of full QCD the chiral behavior, which is an infrared property of QCD, 
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can be severely altered by working on a finite volume. 

For an accurate estimate of the effect of a finite lattice volume on the chiral loops 
one should replace the continuum momentum integrals in the chiral expansion by 
finite sums over the allo wed momenta on the lattice - as already done for baryon 
and vector meson mas ses A full presentation of the results of such an approach 
will be given elsewhere!^ For the present we note that the main consequence of a 
finite volume is to effectively impose a lower limit on the pion momentum integral. 
One can estimate the effect very simply by replacing the factor m 2 /(m 2 + /i 2 ) 
in Eq. |gj by (m 2 + £ 2 )/(to 2 + fi 2 ), where C = 2n/Na = 2tt/L is the lowest 
non-zero momentum available to a pion on the lattice, with a the lattice spacing 
and N the number of lattice sites in each spatial dimension (i.e., a lattice of side 
L fm) . The results of such a simple calculation are illustrated in Fig. 0] from which 
we conservatively estimate that one may need a box L ~ 4 fm on a side to see 
substantial chiral curvature in the extrapolation, even in full QCD. 

3. Bjorken-cc dependence of quark distributions 

In this section we present the results of recent efforts to reconstruct the x depen- 
dence of the quark distributions from their extrapolated moments. We shall briefly 
present updated results for the unpolarized isovector distributions in the nucleon, 
and then focus on new results for the longitudinally polarized valence quark distri- 
bution. An analogous investigation of the valence quark distribution in the pion was 
performed recentlypl whilst analysis of transversely polarized distributions awaits 
the calculation of additional moments. 

3.1. Reconstruction method 

Although the Mellin transforms, Eqs. (J @. have mathematically well-defined 
inverses, such inversion requires knowledge of the behavior of the moments along a 
contour in the complex-n plane. Clearly, evaluation of nucleon matrix elements of 
twist-2 operators (on the lattice or elsewhere) cannot provide this — through the 
OPE, appropriate matrix elements determine moments only for real, integral n. In 
order to proceed, one must assume a parametric structure of the underlying parton 
distribution. This, however, is not a new obstacle; the same problem is encountered 
in parameterizations of PDFs from experimental data at different scales and quark 
distribution functions f{x) (polarized and unpolarized) are often parameterized 
with the form 

xf(x) = Ax b {l-x) c {l + ey/x + -/x). (10) 

The parameters b and c determine the small- and large- a; behavior, respectively, 
and have physical interpretations in terms of Regge behavior and counting rules, 
whereas e and 7 are introduced to provide additional freedom. One easily sees that 
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the corresponding Mellin moments are given in terms of the /3-function by 



(x n ) f =A 



B(l + c, b + n) + eB(l + c, i + b + n) + -fB(l + c,l + b- 



(11) 



Given this parametric form, one can now use the various sets of moments extrap- 
olated from the lattice data to fit the parameters in Eq. Since there are 
at most only three nontrivial moments available for each distribution, one must 
first reduce the number of free parameters by fitting e and 7 to an average of the 
available unpolarizecP^ and polarized^ phenomenological parameterizations. This 
leaves only the parameters A, b and c to be determined from the lattice data. If one 
omits e and 7, it is not clear that b and c retain their physical meanings since there is 
no longer enough freedom in the fits for these parameters to be determined solely by 
the small- and large-x regions, respectively. This is evident when one calculates the 
moments of a known distribution and then attempts to use Eq. to reconstruct 
the distribution f(x). By studying the dependence of the reconstructed parameters 
on the number of moments, N, used in the fit it is apparent that in the full fit, 
Eq. the parameters A, b and c are almost independent of N, while e and 7 
show some variation. On the other hand, if e = 7 = 0, the remaining (purportedly 
physical) parameters show significantly more variation. 

As mentioned in Sec. 2 above, the different crossing symmetry properties of the 
even and odd moments (c/. Eqs. C3>~©) mean that the lattice data correspond to 
valence (odd under charge conjugation, C) moments for n-even (n-odd) in the case 
of the unpolarized and transversity (helicity) distributions, and to the total (even 
under C) moments for the remaining cases. Ideally, one would like several moments 
for each of the C-odd and C-even distributions to reconstruct both the valence 
and the total distributions (which would then allow one to determine the sea). In 
practice, not enough moments are currently known from the lattice, and one must 
be content with reconstructing the valence distributions with the help of additional 
phenomenology. Namely, to recover the valence moments, one needs to shift the 
relevant extrapolated lattice moments by the corresponding moments of the sea 
quark distribution. In practice, since sea quarks are concentrated at small 2, it is 
only the lowest moments (n = in the polarized case, n = 1 for unpolarized) that 
shift to any appreciable degree. For the unpolarized case, the non-singlet sea d — u 
is well determined from Drell-YarP^l and deep inelastic dataj^ (x)d-u — 0.008(1). 
On the other hand, the corresponding quantities for the polarized distributions are 
poorly known and estimates need to be taken into account in the analysis. 



3.2. Reconstructed distributions 

The unpolarized valence isovector distribution reconstructed-^ from the lattice mo- 
ments is shown in Fig. El (solid line). To present realistic error bands on the recon- 
structed distribution, a large ensemble (of order 200 elements) of sets of moments 
is generated by randomly varying each moment within the extrapolated lattice 
error bars. The lightly shaded region then corresponds to the envelope of these 
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n r 




Fig. 5. Reconstructed isovector valence quark distribution x(u v — d v ) in the proton at Q 2 = 
4 GeV 2 . The central fit curve (solid line) and erro r band (lightly shaded) are compared with the 
envelope of the phenomenological distribution^^! (darkly shaded) . 



distributions (the humps in the error band arise from taking the extrema of the 
distributions). The parameter values and errors are quoted in Tableland are given 
by the mean and standard deviation over the ensemble. The darkly shaded regi on 
corresponds to the average of (next-to-leading order) PDF parameterization*^ of 
u v (x) — d v (x) at Q 2 — 4 GeV 2 . The agreement between the reconstructed distribu- 
tions and the phenomenological paramctcrizations is clearly excellent. 



Table 2. Fit parameters for the reconstructed distri- 
butions. For each polarized data set, the upper (lower) 
rows show the parameters for fits with b (c) fixed to its 
phenomenological value. In the unpolarized (polarized) 
fits, e = 1.96 and 7 = 29.1 (e = -0.63 and jj = 10.1) 
were set to their phenomenological valuesl u l I 



Data set 


A 


b 


c 


Unpolarized 


0.21(5) 


-0.63(3) 


3.8(6) 


Polarized, set I 


0.47(3) 
1.9(2) 


-0.51 
-0.05(7) 


2.1(2) 
3.69 


Polarized, set II 


0.28(5) 


-0.51 


1.4(2) 




5(2) 


0.8(5) 


3.69 


Polarized, set III 


0.61(4) 


-0.51 


2.5(1) 




1.5(1) 


-0.24(6) 


3.69 


Polarized, set IV 


0.70(5) 


-0.51 


3.0(2) 




1.1(1) 


-0.39(5) 


3.69 
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In contrast to the unpolarized distribution, there are a number of issues that 
complicate the analysis of the polarized distributions. Firstly, from the cross- 
ing symmetry properties of spin-dependent structure functions, the n — mo- 
ments extracted from the lattice correspond to moments of the total distribution, 
Aq(x) + A.q(x), and in order to construct the valence moment one must subtract 
(twice) the polarized sea moment. Unfortunately, the only experimental informa- 
tion on (1)au-a<2 

from the HERMES datcPH does not provide strong constraints. 
Consequently, we shall investigate two scenarios: (1)aq-A(2 — ( as suggested in 
models in which the non-perturbative sea is generated through meson loop J^j) , 
and (I)ag-aJ = 0.2(2) (as in quark models with SU(6) symmetry, when the Pauli 
exclusion principle is applied^). 
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Fig. 6. Reconstructed isovector valence helicity distribution, x(Au v — Ad v ) at Q 2 = 4 GeV 2 . In 
each panel, the solid (dashed) curve and the corresponding error band is a best fit with b (c) in 
Eq. 1101 fixe d, w hile the darkly shaded region corresponds to the envelope of the phcnomcnological 
distributionsl=il The four sets of curves are described in the text. 



Secondly, as discussed above, finite volume effects have been founds to be par- 
ticularly significant in lattice calculations of qa — (l)Au-Ad, and to account for 
these we shall consider the effects of shifting the extrapolated moments towards 
their physical value. Finally, fitting the three parameters A, b and c with only three 

moments is unreliable (though possible), so one must presently set either b or c to 

loi I 

its phenomenological value. 1 ^ 1 

In Fig. we show the reconstructed PDFs for various shifting scenarios com- 
pared with phenomenological parameterizations. Each panel shows the fits and error 
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bands with b fixed (solid curve, medium shading), and with c fixed (long-dashed 
curve, light shading). The darkly shaded region corresponds to the envelope of the 
phenomenological distributions!^ In the upper left panel, the unmodified extrap- 
olated lattice moments are used (data set I). Fits to these with either b or c fixed 
give a somewhat different reconstructed distribution, and the PDFs extracted from 
the lattice are not in agreement with the phenomenological fits. 

In the upper right plot we show the resulting fits that arise from shifting the 
n = moment up by 0.2(2) to correct for the sea contribution which enters the 
lattice extrapolated moment (data set II) . The agreement between the two fits (with 
b and c held fixed) is still poor and the errors are considerably increased. In the lower 
panels, we contemplate shifts of the moments to correct for possible finite volume 
(and other) effects in the lattice data and their extrapolations. On the left, we shift 
just the n = moment, gA, so that it agrees with experiment (data set III), while 
on the right all three moments are shifted by the same relative percentage towards 
the experimental moments (data set IV). The latter is the best case scenario, with 
the agreement between the two reconstructed distributions significantly improved, 
and (not surprisingly) in excellent agreement with the phenomenological fits. The 
parameters for the various fits are given in Tabled From the variations in Fig. El one 
can conclude that the lattice data currently limit the accuracy of the reconstructed 
x dependence of the isovector valence helicity distribution, and that it is vital to 
understand the current discrepancy between the lattice and experimental values of 
gA- We shall discuss the necessary improvements in Sec. 4. 

3.3. Quark mass dependence of x distributions 

In the moment analyses in Sec. 2 the most dramatic variation with of the mo- 
ments occurs at small values, as one approaches the chiral limit. The transition 
to the heavy quark limit, on the other hand, appears to have little visible effect on 
the moments, and in practice one obtains essentially identical fits to the lattice data 
whether one explicitly imposes the heavy quark limit or not. 

A more graphic illustration of the change of the structure of the nucleon from the 
physical region to the heavy quark limit is provided by the quark mass dependence 
of the x distributions. The chiral extrapolation formulas of Sec. provide moments 
for any value of m w , which allows one to trace how the x dependence changes as 
one goes from the chiral limit to the heavy quark limit, where the distribution 
approaches a <5-function. 

As an example, we consider the unpolarized valence u v — d v distribution, which 
we plot in Fig. for various pion masses from = to m„ = 5 GeV. As m w 
increases, the distribution becomes more sharply peaked, with the peak moving 
towards the limiting value of x = 1/3. This corresponds to fit parameters b and 
c becoming larger, as the x dependence of the valence distribution becomes less 
singular at small x. The distribution at m v = 5 GeV is already rather steep, and 
resembles a constituent quark-like distribution, peaking near x ~ 1/3. Interestingly, 



14 W. Detmold, W. Melnitchouk & A. W. Thomas 




Fig. 7. Dependence of the isovector unpolarized valence quark distribution on the pion mass. 

the shape of this distribution is very similar to that extracted from the moments 
extrapolated linearly to the physical pion massP^ A similar analysis can in principle 
also be performed for polarized PDFs, although the present uncertainties in the m w 
dependence of the polarized lattice moments and the polarized sea prohibit a reliable 
reconstruction of the x dependence as a function of m w . 

4. Conclusions 

The study of the PDFs of the nucleon using lattice QCD has already achieved sev- 
eral notable successes. In the unpolarized case the level of agreement between the 
first three non-trivial moments of the isovector distribution and experiment is ex- 
cellent. To realize this level of agreement one must incorporate the correct leading 
non-analytic behavior of chiral perturbation theory into the extrapolation of the mo- 
ments from the relatively large masses accessible in current lattice QCD simulations 
to the physical regime. Given the normalization condition and an additional three 
moments one can reconstruct the x dependence of the isovector valence PDF quite 
reliably and, as shown in Fig.0 the agreement with experiment is quite impressive. 

The analysis presented here for the polarized isovector quark distribution high- 
lights a number of successes and also areas which demand further study. Amongst 
the successes we note: 

• There has been a clear demonstration of the rather unusual degree of sen- 
sitivity of the polarized moments to the volume of the lattice. Moreover, 
in the case of gA, the agreement with the experimental value definitely 
improves as the volume increases.^ 

• Studies of the chiral extrapolatio n proc edure have clearly established the 
need to include virtual A isobars P^^ Once the A is included there is a 
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significant reduction in the chiral curvature as uij — > 0, compared with 
that found in the unpolarized case. 

For the unpolarized and especially the polarized PDFs there are a number of 
challenges which must be answered in order to reach the stage at which lattice 
calculations rival experimental determinations in their accuracy. In particular: 

• In order to better constrain the chiral extrapolations it is important to push 
the lattice simulations to lower pion masses. Ideally this should occur in full 
(unquenched) QCD, however, quenched, and especially partially-quenched 
simulations would also provide valuable information to guide the extrapo- 
lation. For the unpolarized isovector case, u — d, one needs to confirm the 
predicted chiral curvature of the moments at low pion mass. As suggested 
by the estimates shown in Fig. 01 this will most likely require high statistics 
simulations on lattices of dimensions L ~ 4 fm at a pion mass of order 
300 MeV or lower. 

• To better constrain the functional form of the x-dependence of the PDFs, 
calculations of several higher moments {e.g. n = 3, 4, 5) are necessary. 
This is particularly relevant in the polarized case. Such calculations would 
necessitate the non-perturbative calculation of operator mixing coefficients. 

• Finite volume effects must be e xplo red more thoroughly and taken into 
account in future extrapolations!^ The effects of a finite lattice spacing 

must also be incorporated. 

• In the isoscalar case, one has to resolve the numerical problem of obtain- 
ing a reliable signal for the contribution from disconnected quark loops. 
Furthermore, we note that separations of the sea quark distributions from 
the total isoscalar PDFs, based on lattice data alone, have not yet been 
attempted. They will also require an accurate knowledge of moments with 
n > 3 and evaluation of disconnected contributions. 

With such a program we could expect significant advances in our understanding 
of the quark structure of the nucleon (and other hadrons) over the next few years. 
In particular, we look forward to determinations of the various sea and transversity 
distributions. The methods reviewed here will also prove useful in extracting gener- 
alized parton distributions from recent lattice calculations^! of non- forward matrix 
elements of the various twist-2 operators. 
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